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Abstract. The aperture mass has been shown in a series 
of recent publications to be a useful quantitative tool for 
weak lensing studies, ranging from cosmic shear to the 
detection of a mass-selected sample of dark matter haloes. 
Quantitative analytical predictions for the aperture mass 
have been based on a number of simplifying assumptions. 
In this paper, we test the reliability of these assumptions 
and the quality of the analytic approximations, using 
ray-tracing simulations through a cosmological density 
field generated by very large N-body simulations. We find 
that those analytic predictions which take into account 
the non-linear evolution of the matter distribution, 
such as the dispersion of the aperture mass and the 
halo abundance, are surprisingly accurately reproduced 
with our numerical results, whereas the predictions for 
the skewness, based on quasi-linear theory, are rather 
imprecise. In particular, we verify numerically that the 
probability distribution of the aperture mass decreases 
exponentially for values much larger than the rms. Given 
the good overall agreement, comparisons between the 
observed distribution of the aperture mass and the 
theoretical values provide a powerful tool for testing 
cosmological models. 

Key words: cosmology: theory - cosmology: dark matter 
- cosmology: gravitational lensing - cosmology: large 
scale structure of the Universe - methods: numerical 



1. Introduction 

The gravitational distortion of light bundles from distant 
sources in the universe provides a unique means to 
investigate (the statistical properties of) the intervening 
mass distribution. Being observable through the image 
distortion of the distant faint blue galaxy population, 
this cosmic shear effect offers the opportunity to study 
statistical properties of the large-scale structure. In 
contrast to almost all other methods for investigating the 
large scale structure (LSS) - with CMB being the only 
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exception - no assumption about the relation between 
dark and luminous mass is required. 

Based on the assumption that the intrinsic orientation 
of the background galaxies is random, a net alignment 
in the observed galaxy images can be attributed to 
the tidal field (shear). Hence, the alignment pattern 
of the galaxy images directly reflects the properties of 
the mass distribution. For example, two-point statistical 
measures of the galaxy ellipticities (understood here and 
in the following as two-component quantities, with an 
amplitude and an orientation) can be expressed directly 
in terms of the power spectrum of the mass distribution, 
convolved with a filter function (Blandford ct al. 1991; 
Miralda-Escude 1991; Kaiser 1992; Jain & Seljak 1997; 
Bcrnardeau et al. 1997; Schneider et al. 1998, hereafter 
SvWJK; Kaiser 1998; van Waerbeke et al. 1999; Jain 
et al. 1999, hereafter JSW; and references therein). 
The power spectrum completely characterises a Gaussian 
random density field, and so the two-point statistics, like 
the two-point correlation function of galaxy ellipticities 
or the rms shear in an aperture, suffices to extract the 
statistical information contained in the distorted galaxy 
images. Whereas the earlier of the aforementioned papers 
concentrated mainly on predictions for the cosmic shear 
based on the linear evolution of the cosmic density field, 
it was pointed out by Jain & Seljak (1997) that even 
on scales as large as one degree, the non-linear evolution 
significantly affects the expected amplitude of the cosmic 
shear. 

The non-linear evolution transforms an initially 
Gaussian field into a non-Gaussian one, and thus the 
cosmic shear on small angular scales is expected to 
display significant non-Gaussian features. As pointed out 
by Bernardeau et al. (1997) and SvWJK, the skewness of 
the resulting cosmic shear field is a sensitive measure of 
the density factor £lo, since in quasi- linear perturbation 
theory the skewness is independent of the normalisation 
of the initial power spectrum. In order to define the 
skewness, mass reconstruction algorithms such as those 
developed for cluster reconstructions (Kaiser & Squires 
1993; Seitz & Schneider 1996, Seitz et al. 1998a; Lombardi 
& Bertin 1998a, b; and references therein) can be employed 
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to reconstruct the projected density field. This density 
field, appropriately spatially filtered, can then be used 
to calculate the skewness. In contrast to the two-point 
statistical measures mentioned above, which are defined 
directly in terms of the observable image ellipticities, this 
measurement of skewness is more indirect. This causes 
estimates of the statistical error from the data itself not 
to be straightforward. On the other hand, the aperture 
mass M ap , introduced as a measure for cosmic shear in 
SvWJK, is a scalar quantity directly defined in terms of 
the image ellipticities, and can thus be easily used for 
defining a skewness, as well as a second order statistics, 
the rms of M ap (0), where 9 is the angular scale of the 
circular aperture (definitions are given in Sect. 2 below). 
In particular, in contrast to the two-point ellipticity 
correlation function and the rms shear in an aperture, 
for which the filter with which the power spectrum of 
the projected density field is measured is broad, the 
corresponding filter function for the rms of M ap {9) is very 
narrow, and can be approximated very accurately by a 
delta "function", so that M ap (9) directly measures the 
power of the projected density at wavelength £ ss 4.25/0 
(Bartelmann & Schneider 1999). The skewness defined in 
terms of M ap has been considered in SvWJK. 

As in the evolution of the three-dimensional density 
field, where highly non-linear structures like clusters of 
galaxies form, the projected mass density attains strongly 
non-Gaussian features, e.g., the projection of collapsed 
haloes. As a scalar quantity, the aperture mass is ideally 
suited to probe the full probability distribution of the 
projected mass density; in particular, values of M ap (9) 
far out in the non-Gaussian tail signal the presence 
of massive dark matter haloes. Therefore, peaks in the 
distribution of M ap can be used to search for such haloes, 
independent of their luminous properties (Schneider 1996, 
hereafter S96). Indeed, a first application of M ap to large 
scale structure simulations by Reblinsky & Bartelmann 
(1999) revealed that the detection of dark matter haloes 
through the aperture mass is more reliable in terms of 
completeness and spurious detections and suffers less from 
projection effects than optically selected cluster samples. 
Assuming that highly significant peaks are caused by such 
haloes, one can predict their abundance (i.e., number of 
peaks above a certain threshold per unit solid angle) by 
combining the spatial abundance as predicted by Press & 
Schechter (1974) theory with an assumed density profile, 
such as the universal dark matter profile found by Navarro 
et al. (1996, 1997; combined NFW). This idea has been 
put forward by Kruse & Schneider (1999a; hereafter KS1), 
who found that, depending on the cosmological model 
and the redshift distribution of background galaxies, of 
order 10 such haloes per square degree will be detectable 
in deep ground-based optical images, with a signal-to- 
noise ratio larger than 5. Given the rapid evolution of 
wide-field imaging, a mass-selected sample of dark matter 
haloes is now well within reach. Indeed, a first example 



of a shear-detected mass concentration has recently been 
found by Erben et al. (1999). 

All of these predictions on cosmic shear are made 
using simplifying assumptions in order to make analytic 
progress. JSW tested some of these assumptions using 
ray-tracing simulation through a cosmic density field 
generated by very large N-body simulations; similar tests 
have been carried out by van Waerbeke et al. (1999). 
They found that the major approximations made in 
these analytic treatments, namely the so-called "Born- 
approximation" (which projects the density fields along 
"straight lines"), and the neglect of non- linear terms 
in the propagation equations, are very well satisfied. In 
particular, the twist of light bundles, which vanishes 
identically in the usual analytical treatments, is indeed 
very small. The predictions on the projected power 
spectrum using the approximation for the fully non-linear 
power spectrum of Peacock & Dodds (1996) agree very 
well with the numerical results, whereas the predictions 
concerning the skewness are less precise, indicating a 
breakdown of quasi-linear perturbation theory on small 
scales. 

In this paper, we extend the study of JSW to the 
particular application of the aperture mass statistics. We 
use the same simulations as JSW, resulting in a table of 
shear and projected mass density as a function of angular 
position. From the shear, we can simulate observations 
of the aperture mass as a function of position on the 
"data" field, and investigate its statistical properties. In 
particular, we calculate the probability distribution of 
M ap (9), which we find to be highly non-Gaussian, and 
from that we study the dispersion, skewness and kurtosis 
of the distribution. In agreement with JSW, we find that 
the dispersion is accurately predicted by analytic theory, 
whereas the skewness predictions can differ substantially 
from the numerical results. Of particular interest is 
the kurtosis, since it enters the determination of the 
uncertainty of the dispersion measurement due to cosmic 
variance (SvWJK). We find that the kurtosis is a slowly 
decreasing function of angular scale 9, and attains values 
of ~ 3 even on scales as large as 10'; hence, the 
cosmic variance will be the major source of statistical 
error in the measurement of the power spectrum of 
the projected density field from M ap . The analytic 
predictions on the abundance of significant peaks of M ap , 
and thus presumably of dark matter haloes, turns out 
to be remarkably precise, given the strong assumptions 
made. We confirm the high number density of haloes 
detectable with this method. The shape of the probability 
distribution of M ap in the highly non-Gaussian tail, 
predicted by Kruse & Schneider (1999b; hereafter KS2), 
can also be confirmed to be well approximated by an 
exponential function. 



K. Reblinsky et al.: Cosmic shear and halo abundances: analytical versus numerical results 



3 



2. The aperture mass measure M ap 

In this section, we briefly summarise the properties of the 
aperture mass, i.e., its definition, its relation to the shear, 
and its signal-to-noise ratio. For more details, the reader 
is referred to S96 and SvWJK. 



2.1. M ap statistics 

We define the spatially filtered mass inside a circular 
aperture of angular radius 9 around the point £ by 



M ap (C) := / d^ U{\# - CI), 



(1) 



where the continuous weight function U(d) vanishes for 
i? > 6. If C7 (i?) is a compensated filter function, 



di? i9 U(§) = 0, 



(2) 



one can express M ap in terms of the tangential shear 
7t(£; C) at position £ + C relative to C as 



M ap (C) = J d 2 e 7 t(^;C) Q(\i\), 

(Fahlmann et al. 1994; S96), where 
7t(£;C) = -Re( 7 (£ + C)e- 2i ^), 



(3) 



(4) 



and (f> is the polar angle of £. The function Q is related 
to U by 



i) 2 



d#' U(f) - 17(0). 



(5) 



We use the filter function for 1 = 1 from the 
family given in SvWJK: writing U($) = u(d j 8) / d 2 , and 
Q(0) = q(i9/e)/d 2 , we take 
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u(x) = - (l - x 2 ) ( - -ar 



1 



and 
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g(a;) = -^(1-x 2 ), 

7T 

with = = q(x) for a; > 1. 



(7) 



Signal-to-noise ratio 

An estimate of the shear field 7, and thus of the 
aperture mass M ap (i?) through Eq. (||), is provided by the 
distortions of images of faint background galaxies. The 
complex ellipticity of galaxy images is defined in terms 
of second moments of the surface-brightness tensor (e.g., 
Tyson et al. 1990; Kaiser & Squires 1993). Specifically, we 
use here the ellipticity parameter e (Schneider 1995; Seitz 
& Schneider 1997), which is defined such that for sources 
with elliptical isophotes of axis ratio r < 1, the modulus 
of the source ellipticities is given as \e^\ = (1 — r)/(l + r), 



and the phase of the is twice the position angle of the 
major axis. 

The complex image ellipticity e can then be calculated 
in terms of the source ellipticity e( s ) and the reduced 
shear g = 7(1 — k) _1 by the transformation (Seitz & 
Schneider 1997) 



.(b) 
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1 + g*e 



(s) 



(8) 



This relation is valid only for noncritical clusters. For 
critical clusters, it has to be replaced by a different 
transformation. However, as we are mainly interested in 
the weak lensing regime, the above relation is sufficient 
here. 

It has been demonstrated (Schramm & Kayser 1995; 
Seitz & Schneider 1997) that the ellipticity e of a galaxy 
image is an unbiased estimate of the local reduced shear, 
provided that the intrinsic orientations of the sources are 
random. In the case of weak lensing, n <^ 1, one then has 



(e) = 9 ~ 7 



(9) 



by averaging rt|) with the probability distribution of the 
source ellipticities. 

As for the tangential shear component 7t occurring in 
(^J), a similar quantity for the image ellipticities can be 
defined. Consider a galaxy image i at a position i?,; + £ 
relative to the point £ with a complex image ellipticity 
ej. In analogy to (||) the tangential ellipticity C) °f 

this galaxy is then given by 



C) = -Re (e^ + C) e~ 2i4 



(10) 



where <fii is the polar angle of 

We can now estimate the integral (^) by a discrete 
sum over galaxy images, 



1 



(6) M ap (C) = - y>i(0i;C)W 



(11) 



where n is the number density of galaxy images. The 
discrete dispersion Od of the aperture mass M ap (£) is 
found by squaring ( |ll|) and taking the expectation value 
in the absence of lensing, which leads to 



^ = ^2-<nl**l), 



(12) 



where a 2 = (|e' s -'| 2 ). Performing an ensemble average of 
Eq. (na) leads to the continuous dispersion a c 



Ad Q 2 (tf). 



Finally, the signal-to-noise ratio S at position £ is 
_ M ap (C) _ y/2 Ei^i-X) Qfltfil) 



s(C) 



,1/2 



(13) 



(14) 
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In the simulations, we draw the source ellipticities 
from a Gaussian probability distribution, 

Ps(\e {s) \) = ^ T . 1 7 -2^1 CX P (~^?) > ( 15 ) 

7TO-2 [1 - CX P (-CT e )J V °7 / 

where the width of the distribution is chosen as <r e = 0.2. 
Throughout this paper, we assume that the number 
density of the background sources is n = 30arcmin~ 2 . 

2.3. Analytical work done with M ap 

The aperture mass has been considered in the framework 
of blank field surveys in a variety of earlier publications. 
Introduced as a convenient statistics for cosmic shear, 
SvWJK have calculated the rms of M ap as a function of 
angular scale, using the Peacock & Dodds (1996) approxi- 
mation for the non- linear evolution of the power spectrum 
of density fluctuations. Like other two-point statistics, 
the dispersion of M ap is an integral over the power 
spectrum of the projected mass distribution, weighted 
by a filter function. The filter function corresponding to 
(M ap ) is very narrow and can be well approximated by 
a delta function (Bartclmann & Schneider 1999). Hence, 
(Mf p (#)) reproduces the shape of the projected power 
spectrum and, depending on the cosmological model and 
the redshift distribution of the sources, it reveals a broad 
peak at 6 <~ 1'. One convenient property of the aperture 
mass is that the correlation function of M ap of two 
apertures spatially separated by A9 quickly decreases and 
already achieves values of 1CP 2 for A8 ~ 26. This means 
that measurements of M ap from a large consecutive 
area can be considered independent if the apertures are 
densely laid out on this data field; this is in contrast to 
the rms shear in apertures which is strongly correlated, 
and thus must be obtained from widely separated regions 
on the sky. 

Being a scalar quantity, M ap can also be used 
for higher-order statistical measures of the cosmic 
shear. SvWJK calculated the skewness of M ap , using 
Eulerian perturbation theory for the evolution of the 
three-dimensional density contrast S. In agreement with 
Bernardeau et al. (1997) they found that the skewness is 
a sensitive function of the cosmic density factor Qq, and 
is in this approximation independent of the normalisation 
of the power spectrum. 

A measurement of the dispersion of M ap is affected 
by two main sources of statistical error: the intrinsic 
ellipticity distribution of the source galaxies, and cosmic 
variance. To estimate the latter, one needs to know 
the kurtosis of M ap which cannot easily be determined 
analytically. 

Values of M ap much larger than its rms probe the 
highly non-Gaussian regime of the projected density 
field. From its definition, one sees that large values 
of M ap are expected if the aperture is centred on a 
density peak with a size comparable to the filter scale 



0. Therefore, the aperture map can be used to search 
for such density peaks, presumably collapsed dark matter 
haloes, in blank field imaging surveys. In this way it is 
possible to obtain a mass-selected sample of such haloes 
(S96). Simple analytical arguments in S96 suggest that 
dark matter haloes with an approximately isothermal 
profile are detectable with a signal-to-noise ratio larger 
than 5 if their velocity dispersion exceeds ~ 600 km/s, 
assuming a number density of background sources of 
n ~ 30 arcmin -2 . Indeed, this theoretical expectation 
was verified in the lensing investigation of the cluster 
MS1512+36 (Seitz et al. 1998b). This cluster has a 
velocity dispersion of about <~ 600 km/s, as obtained 
from strong lensing modelling and from spectroscopy of 
cluster members, and is detected in the weak lensing 
analysis with very high statistical significance. 

Assuming that the high signal-to-noise peaks of M ap 
are due to collapsed dark matter haloes, one can attempt 
to estimate the abundance of such peaks using analytic 
theory. KS1 have calculated the number density of 
haloes with aperture mass larger than M ap , N{> M ap , 6), 
assuming (1) that dark matter haloes are distributed 
according to Press & Schcchtcr (1974) theory which 
yields the number density of collapsed haloes as a 
function of halo mass and redshift, and (2) that the 
azimuthally-averaged projected density profiles of these 
haloes can be described by the projection of the universal 
halo density profile found in numerical simulations by 
NFW. Depending on the cosmological model and on the 
redshift distribution of the faint galaxies, the number 
density of peaks of M ap with a signal-to-noise ratio larger 
than 5 was estimated to be <; 10 per square degree, 
and the redshift distribution of these haloes is strongly 
dependent on the behaviour of the linear growth factor for 
density perturbations, and thus on fio- This abundance 
is encouraging, since it allows one to obtain samples of 
haloes selected by their mass properties alone (for a first 
example, see Erben et al. 1999). 

Using the same model, KS2 have calculated the 
probability distribution of M ap for values of M ap much 
larger than its rms, assuming that this non-Gaussian 
tail of the probability distribution is dominated by dark 
matter haloes. They found that the distribution is very 
well described by an exponential; i.e., the tail is much 
broader than for a Gaussian. 

All these analytic predictions are based on a number 
of approximations and simplifying assumptions. In Sect. 4 
below we shall compare these analytic results with those 
found in ray-tracing simulations through a cosmological 
mass distribution obtained from very large N-body 
calculations, as described in the next section. 
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3. Generation of shear maps with ray-tracing 
simulations 

Simulated shear maps due to weak lensing by large-scale 
structure are made by performing ray tracing simulations 
through the dark matter distribution produced by N- 
body simulations (JSW). The N-body simulations used 
are a set of adaptive particle-particle/particle-mesh 
(AP 3 M) simulations. The long-range component of the 
gravitational force is computed by solving Poisson's 
equation on a grid. The grid calculation is supplemented 
with a short range correction computed either by a direct 
sum over neighbouring particles, or, in highly clustered 
regions, by combining a calculation on a localised 
refinement mesh with a direct sum over a smaller number 
of much closer neighbours. The parameters used by the 
N-body simulations are given in Table [l]. 

The simulations were run with a parallel adaptive 
AP 3 M code (Couchman et al. 1995; Pearce & Couchman 
1997) kindly made available by the Virgo Supercomputing 
Consortium (e.g. Jenkins et al. 1998). They followed 
256 3 particles using a force law with softening length 
Z so ft ~ 30 /i _1 kpc at z = (the force is ~ 1/2 its 1/r 2 
value at one softening length and is almost exactly 
Newtonian beyond two softening lengths). Z so ft was kept 
constant in physical coordinates over the redshift range 
of interest to us here. The simulations were carried out 
using 128 or 256 processors on CRAY T3D machines 
at the Edinburgh Parallel Computer Centre and at the 
Garching Computer Centre of the Max-Planck Society. 
These simulations have previously been used for studies 
of strong lensing by Bartelmann et al. (1998), for studies 
of dark matter clustering by Jenkins et al. (1998), and 
for studies of the relation between galaxy formation and 
galaxy clustering by Kauffmann et al. (1999a, b), and 
Diaferio et al. (1999). 

The ray tracing simulations of weak lensing from which 
we use the convergence and shear maps were computed 
by JSW. They used a multiple lens-plane calculation 
that implements the discrete recursion relations for the 
position of a given photon and for the Jacobian matrix 
of the lens mapping at this position (Schneider & Weiss 
1988; Schneider et al. 1992; see Seitz et al. 1994 for 
a thorough justification for this approach). Aside from 
the distance factors, the main input into the recursion 
relations is the shear matrix at each lens plane. The ray 
tracing algorithm consists of three parts: constructing 
the dark matter lens planes, computing the shear matrix 
on each plane, and using these to evolve the photon 
trajectory from the observer to the source. The details 
involved at each step are as follows: 

1. The dark matter distribution between source and 
observer is projected onto 20 — 30 equally spaced (in 
comoving distance) lens planes. The particle positions on 
each plane are interpolated onto a grid of size 2048 2 . 
Since the three-dimensional mass distribution is taken 



from a single realisation of the evolution of the LSS, the 
projected mass distributions of consecutive lens planes are 
correlated. In order to dccorrclatc them, the projection 
is carried out along a randomly chosen one of the three 
coordinate axes; in addition, the origin of the coordinate 
system in each lens plane is translated by a random 
vector and the lens plane is rotated by a random angle. In 
this way, the projected mass distributions of consecutive 
lens planes are as independent as possible, given the 
restriction of only a single realisation of the 3-d matter 
distribution. 

2. On each plane, the shear matrix is computed on a 
grid by Fourier transforming the projected density and 
using its Fourier space relation to the shear. The inverse 
Fourier transform is then used to return to real space. 

3. The photons are started on a regular grid on the first 
lens plane. Perturbations along the line of sight distort 
this grid and are computed using the relation between 
deflection angle and projected density. Once we have the 
photon positions, we interpolate the shear matrix onto 
them and solve the recursion relations for the Jacobian of 
the mapping from the n-th lens plane to the first plane. 

4. Solving the recursion relations up to the source plane 
yields the Jacobian matrix at these positions. Note that 
the ray tracing is done backwards from the observer to 
the source, thus ensuring that all the photons reach the 
observer. The first lens plane is the image plane and 
has the unperturbed photon positions. All sources are 
assumed to be at a redshift of z s = 1. 

There are two kinds of resolution limitations in the 
ray-tracing simulations. The first reflects the finite size 
and resolution of our N-body simulations, the second the 
use of finite grids when computing deflection angles and 
shear tensors on the lens planes. At the peak redshift of 
the lensing contribution, both effects give a small scale 
resolution of order 0.2'. However, since the lens efficiency 
is not very sharply peaked, effects at other redshifts 
also enter. Thus depending on the statistical measure 
being used, the small scale resolution lies in the range 
-0.2' -0.4'. 

On large scales the finite box-size of the N-body 
simulations sets the upper limit on the angular scales 
available. The angular size of our simulation box at z — 1 
is about 3°. Thus on scales comparable to 1°, only a few 
sample regions are available, leading to large fluctuations 
across different realisations. We therefore restrict our 
considerations to apertures with radius 9 < 10' using 
one realisation for every cosmological model. For the 
tCDM model, we use ten different realisations of the ray 
tracing simulations (i.e., they differ in the direction of 
projections, the translation and rotation of the projected 
matter distribution in the individual lens planes) to 
estimate the cosmic variance. 
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Table 1. Parameters of the N-body simulations. 



Simulation 


SCDM 


rCDM 


ACDM 


OCDM 


iVpar 


256 3 


256 3 


256 3 


256 3 


Loitlh^ 1 kpc] 


36 


36 


30 


30 


r 


0.5 


0.21 


0.21 


0.21 


iboxf/i" 1 Mpc] 


85 


85 


141 


141 


n 


1.0 


1.0 


0.3 


0.3 


A 


0.0 


0.0 


0.7 


0.0 


H [km/s/Mpc] 


50 


•50 


70 


70 




0.6 


0.6 


0.9 


0.85 


m p lO lo /i- 1 M 


1.0 


1.0 


1.4 


1.4 


field size [°] 


2.7 


2.7 


3.4 


3.9 



4. Application of M ap to simulated shear maps 

For each of the shear maps generated as described 
in the last section, we create a 2-dim. "M ap map" by 
simulating "observations" of M ap as a function of position 
on the 2-dim. shear maps. The probability distribution 
function of M ap (PDF) and some of its moments are 
then calculated for every M ap map and compared to the 
analytical model. 

It is most instructive to consider two different sets of 
simulated maps: in the first, we neglect noise from the 
intrinsic ellipticity distribution of the background sources 
and compute M ap directly from the shear values on the 
grid according to Eq. (|j|). We do this either in the limit 
of weak lensing, i.e. we use (Q) directly, or we replace 
7t in (||) by the reduced shear gt, which is the quantity 
estimated from the observable galaxy ellipticities. 

In the second set of simulations we introduce elliptici- 
ties of background galaxies according to the distribution 
function (|l5|). The ellipticities add noise to M ap . 

The noise-free results are the ones best compared 
to the analytic results, whereas the ones accounting for 
intrinsic ellipticities yield a more realistic description 
of the observational situation. In the following the 
term "without noise" will refer to the first set of M ap 
simulations, while the term "with noise" will be used for 
the second one. 

As an illustrative example, the 2-dimensional distri- 
bution of M ap for a standard CDM (SCDM) and an 
open model (OCDM) is shown in Fig. [|. In both cases 
high peaks in these maps correspond to haloes in the 
intervening matter distribution. It is possible to construct 
a shear-limited sample of haloes from these maps and to 
determine their abundance. 

Comparing the two model universes, we see that the 
M ap maps reflect the different growth of structure in 
different cosmologies. The M ap map of the OCDM model 
is dominated by many isolated peaks which correspond 
to already collapsed dark matter haloes. The level of 





Fig. 1. The 2-dimensional distribution of M ap for a 
standard CDM (SCDM, upper panel) and an open model 
(OCDM, £l m = 0.3, lower panel), with parameters given 
in Table 1. The field size in both panels is 2°. 



background noise coming from matter not yet collapsed is 
considerably smaller than for the SCDM model in which 
the structure forms later. The peaks in the SCDM model 
are less pronounced and isolated than in the open model. 
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4-1. The PDF of M ap and its moments 

Once we have computed the 2-dimensional distribution 
of M ap , it is straightforward to determine the one- 
point probability distribution function (PDF) of M ap 
and its moments. The PDF contains the cosmological 
information. The lower order moments like rms value and 
skewness can be derived analytically under simplifying 
assumptions, but the PDF itself cannot be calculated. 
Therefore, ray tracing simulations provide the only tool 
for testing the precision of the analytical calculations. 

The qualitative features of the PDF for different filter 
scales 9 and for the four different cosmologies (Table [j]) 
can be studied in Fig. ||. The first point to note is that 
the non-Gaussian features, namely the tail of the PDF 
at high M ap values, are less pronounced for larger filter 
scales. This is due to the fact that the smaller filter scales 
are more sensitive to the already collapsed, non-linear 
objects. The second feature to note is the exponential 
decrease of the tail of M ap which was already obtained 
semi-analytically in KS2. We shall discuss this feature in 
more detail later in this section. 

We now turn to the rms value {M^) 1 / 2 of M ap . Fig. || 
compares the analytical rms value of M ap calculated using 
the nonlinear power spectrum of Peacock & Dodds (1996) 
to the rms values computed from the PDFs without noise 
for 7t (left panel) and gt (right panel). The comparison 
of the latter shows that the difference between shear and 
reduced shear is negligible even on filter scales as small 
as 9 ~ 2 arcmin corresponding to the highly nonlinear 
regime of the mass distribution. 

In the left panel of Fig. ||, there is an excellent 
agreement between the analytic predictions and the rms 
values computed from simulations for the SCDM model. 
There is also good agreement for the ACDM and OCDM 
models, especially for the larger apertures. The notable 
exception is the tCDM model, for which the simulations 
for small filters deviate by a larger factor from the 
theoretical predictions. 

When interpreting this difference between analytical 
calculation and simulation in the rCDM model, one has 
to keep in mind that the numerical results of Fig. ^ are 
based on a single realisation. As the cosmic variance is 
relatively large, it is possible that the large deviation 
is due to the special choice of the realisation. This 
interpretation is supported by the fact that the mean 
for the 10 realisations is considerable lower than for the 
single realisation plotted. Furthermore, the field sizes of 
the simulated fields used are too small to represent a 
characteristic region of the universe. 

The next higher moment of the PDF is the skewness, 
which is defined as 



(16) 
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The skewness S3 of the PDF of M ap as defined 



for which we can perform a similar analysis as for the 
rms value of M ap . As pointed out by Bernardeau et 



Fig. 4. 

in ([l6]) as a function of filter scale 9 for the same 
cosmological models as in Fig. |^. The analytical skewness 
(lines) from quasi-linear theory is compared to the 
skewness obtained from the PDF for both, the tangential 
shear -f t and the reduced shear g t . Errors on the rCDM 
model like in Fig. |[ 

al. (1997), van Waerbeke et al. (1999), and JSW, the 
skewness defined in analogy to (|l^) using a top-hat filter 
is a very sensitive probe of the cosmic density parameter 

n . 

The dependence of the skewness on filter scale 9 is 
displayed in Fig. Again, we compare the skewness 
computed from the PDF obtained from the ray tracing 
simulations without noise, both using 7 t and gt, to 
the skewness of M ap obtained using quasi-linear theory 
(SvWJK). The error bars on the skewness for the rCDM 
model for 2, 5, and 10 arcmin are derived from the 10 
different realisations and are centred on their arithmetic 
mean. 

Again, the differences between the skewness obtained 
from simulations with 7 t and g t are small, though 
slightly larger than for the dispersion, owing to the larger 
contribution from high-K regions to the skewness. This 
difference, which is of order a few percent at most, has 
been predicted to be small in the Appendix of SvWJK. 

When comparing the skewness as determined from 
second-order perturbation theory for the density evolution 
to that obtained from simulations (either computed with 
7t or g t ) we see that the former underpredicts the 
skewness by factors of up to 2. This failure of quasi- 
linear theory for the prediction of higher-order moments 
has been demonstrated previously (Jain & Seljak 1997; 
Gaztanaga & Bernardeau 1998). As we only determine 
the skewness on scales below 10 arcmin, we are in a 
regime where the density contrast is non-linear already. 
The skewness as calculated by Hui (1999) using the so- 
called hyper-extended perturbation theory (Scoccimarro 
& Frieman 1999) may provide a more accurate analytical 
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Fig. 2. The normalised PDF of M ap 
for different filter scales 8 and cos- 
mologies: SCDM (solid line), tCDM 
(dotted line), OCDM (dashed line) 
and ACDM (dashed-dotted line). 
The histograms are obtained from 
M ap maps without noise. Note the 
different scales on the horizontal 
axis. 
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Fig. 3. The rms value of M ap 
computed with the filter (^J) versus 
filter scale 9 for different cosmolo- 
gies. Lines refer to analytic values of 
(Mfp) 1 / 2 from SvWJK, while sym- 
bols refer to rms values obtained 
from simulations without noise, us- 
ing 7 t (left panel) and g t (right 
panel). The error bars in the left 
panel are determined from ( pp| ) and 
©. The symbols for OCDM and 
tCDM are slightly offset on the 
#-axis for better display. The error 
bars in the right panel show the 
standard errors from 10 realisations 
of the rCDM model for 2,5, and 
10 arcmin. They are centred on 
the arithmetic mean (not on the 
realisation plotted). 



prediction of S3 than that from second-order perturbation 
theory. 

Another point to note is the increase of the skewness 
towards smaller filter scales. Generally speaking, such a 
behaviour is expected, as the non-linear structure growth 
becomes more and more important for small filter scales. 
This increase is described insufficiently by quasi-linear 
theory: for the two EdS universes and even for the A 
model on large filter scales above 5 arcmin, this increase 
(not the absolute value!) is predicted satisfactorily, but 
the slope for the open model is larger than analytic 
values on all scales displayed. This discrepancy between 



fully non-linear simulations and quasi-linear theory can 
be attributed to the fact that the open model is much 
more dominated by already collapsed, non-linear objects 
than all other models. 

The highest moment we consider explicitly is the 
kurtosis $4 

(Af 4 ) 

s -^--=h4- 3 - (17) 

The kurtosis is not only important by itself, but also 
for the determination of the error of the rms value of 
M ap , as will be discussed. As for the skewness, the 
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Fig. 5. The kurtosis S 4 of the PDF of M ap [Eq. (|1 
as a function of filter scale 9, for the four cosmological 
models. The kurtosis is derived from the PDF for both, 
the tangential shear j t and the reduced shear g t . Errors 
on the tCDM model are like in Fig. |[ No analytic 
estimate of the kurtosis has been calculated. 



kurtosis for the noise-free simulations for both j t and 
gt is plotted, and the scatter for tCDM is determined 
from the 10 realisations. No analytic result for S4 is 
available; however, using third-order perturbation theory, 
Bernardeau (1998) has calculated the kurtosis for a 
top- hat filter. 

We clearly see that the difference between -f t and 
gt becomes important for the kurtosis, at least for the 
smaller filter scales, since it is even more dominated by 
the non-Gaussian tail of the PDF than the skewness. 
The large error bars on the kurtosis are mainly due to 
large cosmic variance in combination with the small fields 
used; thus, the current simulations are unable to provide 
an accurate determination of S4. 

We now turn to the error bars on the rms values of 
M ap in Fig. ||. In the right panel, they were estimated as 
the standard deviation from 10 different realisations for 
the rCDM model. The error bars in the left panel were 
calculated as follows: 



As shown in SvWJK, an unbiased estimator of (M. 



f2 



from a single aperture is given by 



M 



(tt0 



2\2 



N{N - 



N 



Qi Qj ^ti ^tj 



(18) 



where N is the number of galaxies in the aperture, and 
Qi is the value of the weight function Q for the t-th 
galaxy. The dispersion of this estimator is 



a 2 (M)^S i (M 2 ap ) 2 + 



( 6a e 2 

\5y/2N 



V2(M a 2 p ) 



(19) 



where the two terms in parenthesis correspond to the 
noise from the intrinsic ellipticity distribution, and the 



Gaussian cosmic variance, respectively, whereas the term 
involving S4 is the excess cosmic variance due to non- 
Gaussianity. For a collection of Nf independent apertures, 
all containing the same number of galaxy images, an 
unbiased estimator for (M 2 p ) is the mean M. of M over 
these apertures, and the dispersion is 



a{M) 



a{M) 



(20) 



Note that this result does not assume that the density 
field is Gaussian. If one had a collection of iVf fields 
widely separated on the sky, they would be statistically 
independent, so that Nf — N. In the opposite situation 
where a consecutive area on the sky is available, one 
can lay down apertures on that field, but they will not 
be statistically independent. However, as was shown in 
SvWJK, the M ap values of two apertures which touch 
each other (i.e., with separation twice their radii), are 
almost uncorrelated. Whereas the fact that the two 
aperture masses in these two apertures are uncorrelated 
docs not imply that they are independent (which would 
mean that the joint probability distribution for the values 
of M ap would factorize) - as would be the case for 
Gaussian fields - we assume the statistical independence 
for estimating the effective number of fields Nf entering 
(po|). Thus, the error bars in the left panel of Fig. || 
are obtained from (^0|), assuming that the number of 
independent apertures is Nf = [9/(26*)] 2 , where O is the 
side length of the simulated shear field. 

In contrast, the error bars plotted in the right panel 
of Fig. H for the rCDM model at the three different 
filter scales 9 — 2,5, 10 arcmin are based on 10 different 
realisations of the ray-tracing simulations and allow one 
to obtain a rough estimate for the error from cosmic 
variance. Notice that the error bars are centred on the 
arithmetic mean of the 10 realisations and not on the 
plotted results from a single realisation. 

Comparing the size of the error bars in both panels, 
we see that both methods give errors of the same order 
of magnitude even though the errors estimated from 
the 10 realisations are smaller than the errors from the 
estimator of (M ap ). There are two possible reasons for 
this: first, the effective number of independent apertures 
is probably larger than our estimate given above, so that 
the error bars on the left panel in Fig. |^ most likely 
overestimate the true error. Second, in the calculation of 
the error bars in the right panel, it was assumed that the 
10 realisations are independent; but as argued in Sect. 
U it is possible that the realisations are not completely 
independent. This would lead to an underestimation of 
the cosmic variance. From Fig. || these two competing 
effects cannot be quantified. It should be noted that at 
least on the largest scale plotted, the contribution of 
the intrinsic ellipticity distribution to the error ( f20| ) is 
completely negligible compared to the cosmic variance. 
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Fig. 6. The halo number density 
N{> M a p, 9) (thin symbols and 
thin solid curves) and N(> M ap ,> 
0.6', 9) (thick symbols and thick 
dashed curves) computed without 
noise as a function of the filter scale 
9 for four cosmologies as indicated in 
the panels. Symbols denote results 
from the simulations (O from 7 t , A 
from g t ) whereas the two line types 
display the corresponding analytic 
results from KS1. A signal-to-noise 
ratio S > 5 is used as detection 
threshold for the haloes. Error bars 
in the upper right panel display 
standard errors from 10 realisations 
for tCDM at 2, 5, and 10 arcmin 
(errors for j t offset to the left, errors 
for g t offset to the right) 
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Fig. 7. The same as Fig. || but 
with noise from intrinsic elliptici- 
ties of the sources added. For the 
theoretical model this is done by 
convolving the values from Fig. |^ 
with a Gaussian, with the dispersion 
obtained from the intrinsic elliptic- 
ity distribution (see KS1). Values 
from the simulations are denoted 
by diamonds (thin and thick sym- 
bols). Tangential ellipticities now 
are obtained according to Eq. (|To|). 
A signal-to-noise ratio S > 5 is as- 
sumed for the haloes. Error bars 
centred on the halo abundances are 
standard errors for 7 realisations of 
the ellipticity distribution at 2, 5, 
and 10 arcmin. Error bars in the up- 
per right panel for rCDM are from 
10 realisations for 2, 5, and 10 ar- 
cmin [error for N(> M5, 9) offset to 
the right, error for N{> M 5 , > C, t ,9) 
offset to the left]. 



4-2. Halo abundances 



As already indicated in Sect. 2.3, high signal-to-noise 



peaks of M ap can be identified with dark matter haloes, 
rendering the construction of a mass-limited (more 
correctly: shear-limited) sample feasible. Analytically, the 
halo abundances can be modelled using the Press & 
Schechter (1974) prediction for the mass- and redshift- 
dependent halo number density, and the universal density 
profile of NFW, while in the simulated M ap map all 



connected regions above the corresponding threshold 
are counted as haloes. We shall consider haloes with 
signal-to-noise ratio S larger than 5, i.e., a peak in the 
M ap map is counted as a halo if M ap > M5 = 5<r c ((9). 

We consider two differently constructed halo abun- 
dances in the following: The first sample is simply 
Af(> Ms,9), the number density of haloes with an 
aperture mass larger than M5 for a given filter size 9. 
The second sample selects peaks with M ap > M5 within 



K. Reblinsky et al.: Cosmic shear and halo abundances: analytical versus numerical results 



11 



a connected, cross sectional area of n^, where £ t is the 
corresponding cross section radius; the number density of 
such peaks is denoted N(> M 5 ,> Hence, the size 

of the peaks in the second sample exceeds the threshold 
£ t ; these peaks are expected to be more robust with 
respect to noise coming, e.g., from the intrinsic ellipticity 
distribution and measurement errors. We use a fixed 
value of £t = 0.6 arc minutes. 

In Fig. H the number density of the two halo samples 
as determined from the simulations without noise are 
compared to the results from the analytic calculation 
in KS1 over a range of filter scales 2' < 9 < 10'. The 
four panels in Fig. ^ refer to the four cosmological 
models considered. The error bars for the rCDM model 
at 2,5, and 10 arcmin are again obtained from the 10 
different realisations centred on the arithmetic mean of 
the realisations. 

In general, the number counts determined from 
simulations agree astonishingly well with the analytical 
results, considering the simplifying assumptions entering 
the latter. The deviations between simulations and 
analytical calculation for three of the four cosmologies, 
namely SCDM, OCDM, and ACDM, and especially for 
the filter scales above 5 arcmin, are less than 10 %. The 
largest deviation found for these three models is a factor 
of 2, for the ACDM model at smallest filter scale. 

The only notable exception is the rCDM model 
where the deviation remains above 10 % even for the 
largest filter scales (9 = 10 arcmin). This relatively bad 
agreement has already been noticed for the rms value of 
M ap and is probably due to the fact that the realisation 
plotted is not characteristic for the mean properties of 
that model, as also indicated by the fact that the halo 
abundance lies above the mean of all realisations as 
indicated by the error bars. 

The good agreement between analytic estimates 
and numerical results for the halo number density 
are surprising, given that (a) Press-Schechter theory 
does not exactly reproduce the spatial number density of 
haloes when compared to N-body simulations, and (b) the 
universal density profile found by NFW has been obtained 
by spherical averaging, and therefore cannot account 
for the non-axisymmetry of their projected density. 
Furthermore, (c) the haloes found from the simulated 
M ap are expected to be affected by projection effects 
(Reblinsky & Bartelmann 1999) which are completely 
neglected in the analytic estimates. Despite these effects 
which one might suspect to yield significant discrepancies, 
we find that the analytic estimates are very accurate. 

We also investigate the halo abundance in an 
observationally more realistic situation in Fig. [?], including 
the noise from the intrinsic ellipticity distribution of the 
background sources. The plot displays the same quantities 
as Fig. U except for the fact that the halo abundances 
of the two different samples have been determined using 
the tangential ellipticities in the case of the simulations. 



The analytic estimates are obtained as in KS1. For 
all four cosmologies, we determined error bars using 7 
different realisations of the ellipticity distribution of the 
background sources (|T5j). The error bars from the 10 
realisations shown for rCDM are slightly sub-Poissonian, 
as in Fig. ^. As expected from the large value of the 
kurtosis the error coming from the intrinsic ellipticity 
distribution is much smaller than the error coming from 
the cosmic variance. On the whole, the number of 
detected haloes is increased in all cosmologies because, 
due to the steepness of the Press-Schechter mass function 
for massive objects, there are more objects just below the 
threshold than above it. So on average more objects will 
be lifted above the threshold by noise than brought down 
below it. 

4.3. The tail o/M ap 

In Fig. H we compare the PDF for M ap > M5 as obtained 
from analytic calculations (KS2) with that derived from 
the simulations without noise. The PDF is shown for 
four filter scales 9 = 2,4,6,10 arcmin in the range 
M5 < M ap < 2M5, for which the analytic results predict 
a nearly exponential behaviour. Indeed, the numerical 
PDF in the non-Gaussian tail also seems to follow an 
exponential rather closely, with a slope very similar to 
the analytic result. 

In order to see how much the PDF varies between 
different realisations, we have plotted in Fig. ^ the PDF 
for M5 < M ap < 2M5 obtained from the 10 realisations 
in the rCDM model, for 3 filter radii, together with their 
mean and the corresponding analytic prediction. We find 
that for the smallest filter scale 9 — 2' , all 10 realisations 
are clearly below the analytic result, whereas for the 
larger filters, the realisation mean of the PDF agrees very 
well with the analytic prediction. 

Remembering that the analytic predictions were made 
by assuming that all high values of M ap are coming 
from regions close to collapsed haloes, in addition to 
the assumptions used for estimating the number density 
of M ap peaks (Press-Schechter halo abundance and 
NFW density profile), this good agreement is somewhat 
surprising. 

5. Conclusion 

We used ray-tracing simulations through N-body- 
generated cosmic density distributions to study the 
statistical properties of the aperture mass M ap as a 
statistics for cosmic shear measurements and for finding 
dark matter haloes from their shear properties. In par- 
ticular, we have compared results from these simulations 
with the available analytic results and found in most cases 
a very good agreement, except for the skewness which is 
the least accurate of these predictions. Whereas all other 
predictions tested here are based on manifestly non-linear 
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Fig. 8. The tail of the PDF of 
M ap for the same cosmologies as 
indicated in Fig. || for different filter 
scales. In each panel we plot the 
PDF obtained from the analytic 
estimate in KS2 with thin lines 
and that from simulations with 
thick lines. The line types spec- 
ifies the cosmology: SCDM (solid 
line), rCDM (dotted line), OCDM 
(dashed line) and ACDM (dashed- 
dotted line). The M ap -range is 
[M 5 , 2M 5 ], where M 5 = Sx 0.016/6* 
with S > 5 for the parameter speci- 
fied in the text. 



results (like the Press-Schechter halo abundance and the 
Peacock & Dodds power spectrum), the skewness was 
estimated analytically by using second-order Eulerian 
perturbation theory which, on the scales considered, is 
not very accurate. 

Comparing the results from our simulations with 
analytic studies, we obtain the following main results: (1) 
The rms of M ap is accurately described by analytic results 
if the fully non-linear prescription of the power spectrum 
of density fluctuations is used. (2) The statistical error 
of this rms is dominated by cosmic variance, which in 
turn depends on the kurtosis of M ap . This kurtosis turns 
out to be relatively large even on angular scales of 
~ 10', implying the need for many more measurements 
of M ap than expected for a Gaussian field, for a given 
accuracy of the estimated projected power spectrum. (3) 
The skewness is only approximately described by analytic 
considerations based on second-order perturbation theory. 

(4) The predicted abundance of dark matter haloes 
detectable at given statistical significance is very well 
approximated by the semi-analytic theory which combines 
the Press-Schechter number density of haloes with the 
universal density profile of Navarro, Frenk & White. 

(5) Similarly, the functional form of the probability 
distribution of M ap for values much higher than the rms 
(i.e., in the non-Gaussian tail) is found to closely follow 
an exponential form, of similar slope and amplitude as 
predicted by analytic theory which needs to assume that 
such high values originate due to collapsed haloes. Thus, 
on the whole, we find that the analytical estimates for the 
statistical properties of M ap are surprisingly accurate. 

We also find that our simulations are not sufficiently 
large for an accurate estimate of the higher-order 
statistical measures, owing to the finite size of the 



simulation box in combination with the large effect of 
cosmic variance. 

As discussed in SvWJK, KS1, KS2, van Waerbeke 
et al. (1999) and Bartelmann & Schneider (1999), 
the aperture mass is a useful cosmic shear measure 
which will eventually allow one to constrain cosmological 
parameters, completely independent of any assumption 
on the relation between mass and light. For this purpose, 
the predictions from cosmology must be known precisely, 
and our results here indicate that analytic estimates 
are relatively accurate. Unfortunately, we found a large 
cosmic variance; e.g., in the estimate of the variance 
of the rms value of M ap , the kurtosis enters and it 
decreases only rather slowly with increasing filter scale. 
It can be expected that the first successful application 
of the aperture mass will be the definition of a sample 
of haloes defined in terms of their lensing properties 
only, with a first example given by Erben et al. (1999). 
The combination of cosmic shear information and CMB 
measurements can be extremely useful, as shown by 
Hu & Tegmark (1999), increasing the precision of the 
determination of cosmological parameters substantially 
over each of the two individual methods. Their study 
was based solely on the dispersion of cosmic shear, i.e., 
on second-order statistics. It is to be expected that a 
similar combination of CMB results with the PDF of 
M ap will yield even more precise parameter estimates. 
A detailed study of this combination is expected to be 
very valuable, but requires a larger grid of cosmological 
N-body simulations. 
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